Averaged dynamics of two-phase media in a vibration field 
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The averaged dynamics of various two-phase systems in a high-frequency vibration field is studied 
theoretically. The continuum approach is applied to describe such systems as solid particle suspen- 
sions, emulsions, bubbly fluids, when the volume concentration of the disperse phase is small and 
gravity is insignificant. The dynamics of the disperse system is considered by means of the method 
of averaging, when the fast pulsation and slow averaged motion can be treated separately. Two 
averaged models for both nondeformable and deformable particles, when the compressibility of the 
disperse phase becomes important, are obtained. A criterion when the compressibility of bubbles 
cannot be neglected is figured out. For both cases the developed models are applied to study the 
averaged dynamics of the disperse media in an infinite plane layer under the action of transversal 
vibration. 

PACS numbers: 47.55.Kf, 47.55.dd, 43.40.At, 47.35.-i 



I. INTRODUCTION 

Oscillatory processes can be found in a sheer uncount- 
able number of situations in nature at various time and 
length scales: from subatomic to astronomic scales. Vi- 
bration is a mechanical oscillatory process with an am- 
plitude small compared to the characteristic length scale 
of the system. Often, as in our study, it is assumed that 
the inherent time scale of the system is much larger than 
the period of the oscillation. Vibration mechanics has 
been studied for a long time, general concepts of this 
nonlinear phenomenon have been put into practice. It 
is widely used in industry for transportation and separa- 
tion of granular matter such as powders and grains, to 
enhance mixing and reaction rates of different species. 
In medicine, ultrasound is applied for diagnostics and 
healing. In microfluidics and biotechnology it serves as 
a means to manipulate physical (colloidal particles, fine- 
dispersed powders, liquid drops in microemulsions) and 
biological (cells, microorganisms, and micromolecules) 
objects suspended in liquids. Quite a broad spectrum 
of applications has motivated a great interest in the fun- 
damental problems of vibration dynamics. The present 
work is aimed at studying the dynamics of two-phase me- 
dia in a high-frequency vibration field. 

There has been much effort to understand the impact 
of vibration action on the behavior of systems of differ- 
ent nature. Particularly, it has been understood that 
vibration can substantially influence the dynamic and 
rheologic properties of systems. An outstanding exam- 
ple is the Kapitza pendulum with a periodically mov- 
ing point of support, "'^1 where the upper vertical state 
of the pendulum becomes stable. Other known effects 
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are a metal ball that lifts in ambient sand medium 
under vibration,^ granular matter that liquidizes'^ or 
demonstrates the formation of oscillons^i in a vibration 
field. As typical hydrodynamic examples we mention 
excitation of a steady relief at the oscillating interface 
pure liquid-particle suspension, ^'^ and thermal vibration 
convection,^ where the action of vibration on the nonuni- 
form fluid results in generation of a slow macroscopic 
motion even in the absence of gravity. 

The behavior of a single solid inclusion in fluid environ- 
ment under vibration has attracted much attention since 
long time ago. As far back as 1831, Faraday used par- 
ticles to visualize standing acoustic waves and observed 
the effect of particle localization in the acoustic fieldi^ 
The dynamics of a solid particle suspended in a verti- 
cally oscillating fluid was studied in Ref. 0. It was shown 
that vibration results in an average force, which can be 
used to affect particle dynamics. Granat-ft considered 
theoretically a motion of a solid sphere in a periodically 
oscillating flow of a viscous fluid. For the fluid with the 
density higher than that of the sphere, the latter oscil- 
lates with a smaller amplitude and a lag in phase. Con- 
versely, in a relatively less dense fluid the oscillation of 
the sphere occurs with a larger amplitude with respect 
to the flow and leads in phase. In an experimental study 
by Chelomey,^^ the behavior of solid bodies in a verti- 
cally vibrated vessel with a fluid was investigated. It was 
observed that under vibration action, bodies with higher 
than ambient density could rise to the surface, whereas 
light bodies could sink. Theoretically, the averaged be- 
havior of a solid sphere suspended in oscillated fluid near 
a wall was discussed in Ref. |3 and has been recently 
generalized. It has been demonstrated that in addi- 
tion to the gravity and Archemedian forces there appears 
an averaged force that attracts the body to the wall due 
to vibration. This vibration force can lead to ascent of 
heavy or settling down of light bodies, as was observed 
experimentallyiii 
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The dynamics of a single compressible bubble in vibra- 
tion fields has been studied by many authors. The first 
analysis of the bubble dynamics was made by Rayleighi^ 
(see also Ref. 15). A nonlinear equation that couples 
the radius of a bubble to the pressure in the liquid at 
large distance from the bubble, known as the Rayleigh- 
Plesset or Rayleigh-Lamb equation, was obtained and 
the eigenfrequencies of the small oscillations were de- 
termined. Since then, many important aspects of non- 
linear bubble dynamics and cavitation hasve been stud- 
ied, e.g., natural and forced oscillations, parametric in- 
stability, which causes the oscillation of the bubble shape 
against a background of spherical pulsation, the infiuencc 
of heat exchange between bubbles and ambient fiuid and 
nonideality of the gas filling the bubble, the behavior of 
bubbles in nonuniform flows (see reviews in Refs.^33)- 

In the context of the time-averaged dynamics of a bub- 
ble an important problem is obtaining the force exerted 
on the bubble due to vibration. The behavior of a single 
bubble in a standing acoustic wave is determined by "the 
primary Bjerknes force" after Bjerknesir^ who proposed 
its qualitative explanation. The effect is quite general: 
the bubbles migrate to the nodes of the pressure wave at 
high frequencies w > and to antinodes at low frequen- 
cies w < where luq is the frequency of natural oscilla- 
tion of the bubble. The average force on a particle, called 
the force of radiation pressure, was obtained by Yosioka 
and Kawasimaii^ The considered setup comprised a liq- 
uid drop in another liquid in a standing acoustic wave, 
where both liquids were assumed to be compressible. A 
generalization for the case of an arbitrary nonuniform 
oscillating flow has been provided by Alekseevf22, These 
results can be applied to analyze both limiting cases - 
the primary Bjerknes effect for deformable bubbles and 
the behavior of a solid particle in nonuniform pulsating 
flows. 

Generalization of the description of a separate particle 
suspended in a fluid to the case of a disperse medium (en- 
semble of many particles) requires space averaging. The 
precise description of all the disperse particles becomes 
redundant and unfeasible. There are basically three pos- 
sibilities to obtain a reasonable model: (i) statistical ap- 
proach, concerned with averaging over the ensemble of 
particles; for solid particles and bubbles such a proce- 
dure is carried out in Refs. I2lll22l2^ (ii) derivation of 
a kinetic equation (see, e.g., Refs. 24*25*26*); (iii) contin- 
uum models, where after the space-averaging the disperse 
particles are treated as a separate continuum. Sij^S In our 
study, we apply the last approach. 

Although, wave propagation in disperse media is well 
understood, the problem of the time-averaged dynam- 
ics for such systems has received relatively little atten- 
tion and remains rather unexplored. Most of studies 
deal with acoustic vibration, when the compressibility 
of carrier fluid is of crucial importance; the feedback of 
particles on the flow is conventionally neglected. For in- 
stance, Ganiev and Lapchinsky^^ studied averaged col- 
lective behavior of bubbles and solid particles suspended 



in fluids in acoustic fields. The analog of the primary 
Bjerknes effect was observed and simple explanation was 
provided. In the present study, we concentrate on the 
averaged effects in various disperse media under the ac- 
tion of high frequency, but subacoustic vibration field, 
i.e., when the fluid is incompressible. We also take into 
account the feedback of particles on the carrier fluid and 
demonstrate that this is an important factor, which can 
lead to nontrivial effects. 

The paper is outlined as follows. In Sec. ^1 the the- 
oretical model is developed for a suspension of nonde- 
formable particles. The equations describing the pulsa- 
tion and averaged dynamics are obtained. On the basis 
of this model the stability of a quasiequilibrium state in 
a plane layer is analyzed in Sec. lIIII Section llVl addresses 
the vibrational dynamics of bubbly fluids. An averaged 
model accounting for the compressibility of bubbles is ob- 
tained. Particularly, it is shown that in the limit of weak 
compressibility this model reduces to the model for non- 
deformable particles, obtained in Sec. ^ In Sec. the 
developed model is applied to investigate the dynamics 
of a bubbly fluid in a plane layer under vibration action. 
The results are summarized in Sec. I VII 



II. THEORETICAL MODEL OF A 
NONDEFORMABLE PARTICLE SUSPENSION 
IN A VIBRATION FIELD 

A. Governing equations and basic assumptions 

Consider the behavior of an isothermal fluid (liquid or 
gas) laden with disperse particles (solid particles, small 
droplets of another liquid or bubbles) under the action 
of high frequency vibration. We assume that all the 
particles are monodisperse spheres of a radius and 
the amplitude of vibration a is small compared to this 
size: a <^ (for viscous pulsation a milder restriction 
is enough to be held, see Sec. Ill B|l . The frequency of 
vibration uj is assumed to be so high that the size of a 
viscous boundary layer £ = y/ v/ui [v is the kinematic 
viscosity of the fluid) near the walls of a container with 
the suspension is small with respect to the characteristic 
scale L of a flow. On the other hand, the scale L is con- 
sidered to be small compared to the acoustic wave length 
Aa — 2'!TVa/oj (va is the acoustic speed), which leads to a 
twofold inequality i//L^ <^ lo <^ Va/ L. 

We start with a two-fluid model valid for dilute sus- 
pensions, when one can neglect interparticle interactions 
and interactions between the particles and walls of a 
container. We assume that the size is large enough 
to neglect diffusion of particles, but small with respect 
to L. On a scale much greater than the particles 
are regarded as a continuum with the volume fraction 
If — A/3TTr^n (the volume fraction of the fluid phase is 
l~ip), where n is the number of particles per unit volume 
of the medium. It is also supposed that the particles can 
neither deform nor combine into agglomerates; the car- 
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rier and the disperse phases have constant densities p and 
Pd (hereafter, the subscript "d" stands for the disperse 
phase). Since if differs from n by a constant factor, in 
the following ip is called concentration. After averaging 
over space the equations for mass and momentum of the 
phases readt^ 

div {l-ip)V + div ipUd = 0, (la) 

^ + </7divUd - 0, (lb) 

(l-(^)p^ = -VP + mV- (l-(^)£-(pF, (2a) 



d\5d 

dt 



(2b) 



where U and Ud are the velocities of the phases, P is 
the pressure in fluid, p, is the fluid viscosity, the shear 
rate tensor is given by e.y — ViUj + S/jUi. Derivatives 
D/Dt = d/dt + V-V and d/dt = d/dt + Vd-V are taken 
along the Largangian path of the elements of the fluid 
and the disperse phase, respectively. 

The interphase interaction force can be written as 

f0ll0WS:28'3"-31 




1 dW 



'2^ dt 



(3) 



Here, the relative velocity of phases W = — U and 
the ratio of viscosities of phases k = Pd/p are introduced. 
The first term in @ is a force due to pressure gradient, 
the next forces are the Stokes (Hadamard-Rybczynski) 
drag^ and three unsteady forces: the Basset history 
force, the added mass force, allowing for the inertia of 
the surrounding fluid, and a second memory force that 
is essential for non-solid particlesi^ We do not consider 
the effects that can be initiated by rotation of particles, 
therefore the Magnus force is not taken into account. 

The relative contribution of the Stokes drag Fs, the 
Basset Fg, and the added mass Fm forces to the inter- 
phase force is governed by an "unsteadiness parame- 
ter" K — rdl\fvtl'- 

where is the timescale of the variation of the relative 
velocity of phases. Note, that the pressure gradient term 
and the added mass force in (O are of the same order. 
The same argument can be applied to the memory forces, 
except for the limiting cases of an inviscid bubble (k = 0) 
and a solid particle (k ^ 1), when, either the Basset term 
or the second memory force vanish, respectively. 

For quasisteady flows (K <C 1) and the flow regimes 
with finite if, the result |j3Jl holds provided that the par- 
ticle Reynolds number Re^ is small and gradients of the 



velocity U are not too largei^ 

TdW 



Red 



< 1, 



vL 



< 1. 



(4) 



In the case of substantially unsteady flows {K ^ 1), the 
force © has inviscid nature and the validity criterion 
does not depend on the viscosity. The restrictions Q 
should be replaced with the conditions 



Wt 



rd L 



(5) 



which arise from linearization requirements. Note, that 
in this limit. Red is not necessarily small and for finite 
and large Red the forces depending on viscosity in ex- 
pression lO have to be modified. However, for the flow 
with X 3> 1, these terms (irrespective of the modifica- 
tion) are small compared to the leading pressure gradient 
and added mass forces and therefore can be neglected. 

The governing equations are formulated in an iner- 
tial reference frame; the infiuence of vibration should be 
taken into account through boundary conditions. 



B. Pulsation velocities 

The presence of two considerably different time scales, 
namely, the characteristic oscillation time 1/uj and the 
dissipative time scale L'^ /v, makes reasonable to use the 
multiscale technique*^ the "fast" and the "slow" aver- 
aged motions are treated separately. So, the operator of 
time differentiation is replaced by a sum of two operators 
d/dt — ^ uid/dr + d/dt, where t = ojtis the fast time; the 
actual fields of the velocities, the pressure, and the con- 
centration are represented by sums of the averaged (over 
the fast time r) and the fast, or pulsation, components, 
respectively: U = u + V, Ud = Ud + Vd, P = p + q, 
ip = c + (p. 

Let us substitute this ansatz into the governing equa- 
tions and obtain the pulsation parts of the momentum 
Eqs. (O. We retain only the leading terms and take 
into account a smallness of the concentration pulsation, 
(j> ^ {a/ L)c, which follows from Eq. (jl b|) . The dimension- 
less variables are introduced using the following units: 
l/u! and L^/iy for the fast and the averaged time, auj 
and v/L for the pulsation and the averaged velocities 
of the phases, pauP'L and pv^ jl? for the pulsation and 
the averaged pressure, respectively. It is assumed that 
the space derivatives can be evaluated with the scale L, 
which means that the small scale dynamics in boundary 
layers near solid walls and turbulent pulsations are not 
considered. As a result, we arrive at the equations for 
the fast dynamics: 



div(l -c)V + div cVd = 0, 
^-HdivcVd = 0, 



(6a) 
(6b) 
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dV 

1 /2/3- 



-V^ + cf, 
9t ' 



(7a) 
(7b) 



K"^ \ l + K 
1 



-i— + f 



gw(C) rfC 



(8) 



where the relative pulsation velocity w = — V and 
a dimensionless parameter 5 — Pdl P-, which is the ratio 
of densities of the phases, are introduced; the pulsation 
component of the memory force F„i is denoted by f^- For 
the case of pulsation motion, the unsteadyness parameter 
^ = Tdjt has the meaning of the ratio of the particle 
radius to the size of a viscous boundary layer around it. 

The restrictions for the amplitude of vibration a are 
imposed by the inequahties ©, and the necessity of 
linearization of the pulsation problem ©-((Hll- For the 
motion with finite K and in the limit of inviscid pul- 
sations, X S> 1, the vibration amplitude is required to 
be small in the sense a <C r^. However, in the limit 
of viscous pulsation, if <C 1, a much milder restric- 
tion for the vibration amplitude is enough to be held: 
a <C min (L, K~^rd) ■ 

Let us restrict our consideration to the case of the 
monochromatic vibration and proceed to complex am- 



plitudes: J^{r, t) = Re ^(r) exp(2r) 



For such a vi- 



bration process, f coincides (in dimensionless form) with 
the Fourier-transform of the pu lsation force f and can be 
represented as (cf. Refs. l3ll33j) : 



2K^ 



1 + A + 



(1 + A)V(A) 



K [A3 - A2 tanh A - 2/(A)] + (A + 3)/(A) 



(9) 



,(10) 



where A = (1 + i)K/V2 and /(A) = A^ tanh A - 3A + 
3 tanh A. 

As follows from Eqs. ©, Q and relation ©, the 
pulsation velocities are not independent, therefore the 
problem can be expressed in terms of an auxiliary field 
V = (1 — c)V + cVrf having the meaning of the amplitude 
of the volume-weighted velocity of the medium: 



VxS'v^O, divv 
XPni + iS 



0, 



s 



(11) 

(12) 



where p„i = 1 — c + cS is the density of the medium and 
l/pe = 1 — c + c/S is the specific volume. The pulsation 



velocities expressed in terms of v are as follows: 



V 



Vd - 



X + iS/pe 

X + i 
X + "iS/pe 



(13) 



(14) 



Note, that except for the limiting cases of viscous {K ^ 
1) and inviscid {K ^ 1) vibration regime, there appears 
a phase shift in the pulsation velocities of the particles 
and carrier fluid. 

Now restrict our theory to the case of a closed non- 
deformable cavity, filled entirely with the two-phase 
medium. Equations Hll|l should be supplemented with 
the boundary conditions. The smallness of the vibration 
amplitude allows us to prescribe the boundary conditions 
at the average position. Because of incompressibility of 
the medium, the normal component of the velocity v has 
to coincide with that of the boundary velocity ^ at every 
point of the boundary T: 



n 



v|r = n • I, 



(15) 



where n is the unit normal to the average position of the 
boundary. 

Assume that the system is subjected to translational, 
linearly polarized vibration. In this case the boundary 
velocity of the cavity is a constant vector: ^ = j (j is the 
unit vector aligned along the vibration axis). 

Since the particle concentration is small, the solution 
to the pulsation problem (|ll|l . H15() is sought as a series 
in c: V = vq + vi + . . . with vi ~ c (and analogously for 
V, Vd). Then, to the leading order we obtain 



vo = Vo = j. 



(16) 



i.e., the pulsation field is constant in the entire cavity: 
uniform fluid moves together with the cavity as a solid. 
Hence, to this approximation, the vibration does not in- 
fluence the system, therefore it is necessary to obtain a 
nonuniform correction to the pulsation field. 

Taking into account Eqs. Hill) . (|15|l and the result l|16|l . 
to the next order we arrive at a problem 



V X vi 
n • vi|r 

Si 



Sijx Vc, divvi = 0, 
0, 

X + i 



(17) 



= (5-1) 



X + iS' 



where 5*1 is a coefRcient in a corresponding expansion of 
(inj): S ^ l+Sic+Oic^). It becomes clear from Eq. l(T7|l . 
that a nonuniformity of the pulsation field vi is caused 
by interaction of vibration with inhomogenieties in the 
particle concentration. We make a note, that an explicit 
solution for the field vi is not needed while obtaining 
the averaged equations of motion: only the vorticity and 
divergency of vi are further used. 
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C. Averaged equations. Single-fluid approximation 

We now turn to averaging of the governing equations. 
Let us write out the averaged parts of Eqs. (Q), Q and 
retain only the leading terms. Note, that averaged dy- 
namics can be described in terms of a single-fluid model. 
Contrastingly to the pulsation motion, relative contri- 
butions of the Stokes drag, the memory forces, and the 
added mass force to the averaged dynamics are defined 
by a small parameter r^/L; the unsteady forces turn out 
to be small compared with the Stokes force and there- 
fore can be neglected. Because the difference in the av- 
eraged velocities of the phases — u is of order r^/L^, 
we neglect it everywhere, except for the interphase force. 
In the latter, this difference must be retained: despite 
its smallness, it is multiplied by a vibration parameter 
Rv = a^uj'^L'^ /v'^ , which can be large, and a correspond- 
ing contribution in the averaged equations is generally 
finite. As a result, we obtain the following equation of 
motion 



dn 



( — -f u • Vu ) + i?,V • [(1 - c) VV + &cWW^ 



= -Vp + V • (1 - c)e. 



(18) 



Here, e is the averaged shear rate tensor; overlines are 
used to denote averaging over time. The accepted ap- 
proximation means, that the particles are "frozen" in the 
averaged carrier flow. 

It is important to note, that the nonuniformity of the 
pulsation field can be caused not only by the considered 
mechanism, when vibration interacts with nonuniformity 
of the medium due to admixture. Another reason for a 
nonuniformity in the pulsation field is interaction of vi- 
bration with the vorticity of the averaged flow. A contri- 
bution to the averaged equations due to the second mech- 
anism is defined by the vector of pulsation transport (see 
discussion in Ref. i7|). However, for the linearly polarized 
vibration, the case under consideration, this vector van- 
ishes. A quite similar situation takes place in thermal 
vibrational convection.^ The only difference is that there 
the nonuniformity of density is caused by nonisothermal- 
ity of fluid, but not due to the presence of a disperse 
phase. 

Using (|13|l . (|14|l . we calculate the averaged Reynolds 
stress 



(1 - c) VV + (JcVrfVrf G Re (vv* 





" 1. 


X + « 









(19) 
'.(20) 



Taking into account smallness of the particle concentra- 
tion, the result (|16|) . and Eqs. dJ, we evaluate the 
"vibration force" 



where dots denote the omitted gradient terms. Such 
terms lead only to a renormalization of the pressure, 
which is not important in the framework of our consid- 
eration. 

In the vein of the Boussinesq approximation, we ne- 
glect the nonuniformities of the density everywhere ex- 
cept for the the vibration force, where the small con- 
centration is multiplied by the vibration parameter i?^,, 
which is generally large. Such a transformation corre- 
sponds to a limiting transition, when i?„ — > c», c — *■ 0, 
but their product i?^c* is finite (here is the reference 
value of the particle concentration). 

Averaging of the mass balance equations is straight- 
forward. Finally, we end up with a set of the averaged 
equations describing the behavior of a suspension of non- 
deformable particles in a high frequency vibration field: 

^ + u • Vu - i?J j • Vc = -VH + V^u, (22) 
at 

dc 

— +\i-Vc = 0, divu = 0, (23) 
at 

where H is the renormalized pressure, the concentration 
is normalized by c* , and the vibration parameter is rede- 
fined Ry = a?u}'^L'^Bc^/v'^ . 

The vibration force, appeared as a result of averag- 
ing of the equations, is caused by interaction of the 
monochromatic linearly polarized vibration with the in- 
homogeneity of medium due to a nonuniform distribu- 
tion of particles; for the uniform distribution of particles 
it vanishes. The force is directed along the vibration 
axis and takes on maximal values in the domains of the 
flow where the gradient of particle concentration coin- 
cides with the direction of the vibration axis. For the 
initial state with a nonuniform distribution of particles 
and the vanishing averaged velocity the action of the vi- 
bration force can result in generation of averaged motion. 

It is important to make a note of the generic character 
of the vibration force. The dependence on the parame- 
ters (5, K, A', i.e., on the the nature of the particles, is 
gathered in the parameter B. This allows us to make 
two conclusions. First, since B > Q for all values of these 
parameters, the direction of the force is determined only 
by the concentration field and an orientation of the vi- 
bration axis. Second, for different media only the inten- 
sity of response towards vibration action changes. How- 
ever, such a quantitative distinction is worth discussing 
for comparison with experiments, therefore properties of 
the function B = B{S, n, K) are discussed in Appendix. 

The set of averaged Eqs. (|23l, must be supple- 
mented with boundary conditions for the velocity. For 
the considered case of vibration action, the Schlichting 
mechanism.'^^ of surface generation of averaged flow is 
absent, therefore on a rigid surface F one should use the 
conventional no-slip boundary condition 



R„V ■ [GRe(vv*) 



B = 



{6 — if' Imx 
'2\x + ^5\^ ' 



(21) 



u r 



0. 



(24) 



Indeed, the intensity of the averaged flow caused by 
the Schlichting mechanism is governed by the pulsation 
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Reynolds number^ Rep = a^oj/v. Compared to the con- 
sidered bulk mechanism, which is governed by the pa- 
rameter Ry, the Schlichting mechanism can give a finite 
contribution to generation of averaged motion only if the 
leading part of the pulsation velocity is nonuniform. In 
our case, this requirement is not satisfied: according to 
(|16|l . the leading part of the pulsation velocity is uni- 
form, a nonuniformity occurs only as a small correction 
to the uniform term. Hence, the Schlichting generation of 
averaged fiow is negligible compared to the contribution 
of the bulk mechanism. Rather similar situation takes 
place in thermal vibration convection: for the same kind 
of vibration action on weakly nonuniform fluid, when the 
density nonuniformity is caused by nonisothermality of 
fluid, the Schlichting mechanism does not contribute to 
generation of average flow. 

III. STABILITY OF A QUASIEQUILIBRIUM IN 
VIBRATION FIELD 

A. Basic quasiequilibrium state 

Let a suspension fill an infinite plane layer — 1 < z < 1 
(of a dimensional width 2h) in weightlessness. Consider 
the infiuence of the transversal linearly polarized vibra- 
tion on the stability of a quasiequilibrium. A quasiequi- 
librium is thought of as a state in which the averaged 
velocity vanishes in the entire volume of the system, 
whereas the pulsation velocity is different from zero. 

This problem is described by Eqs. it^ . jS^J with j = 
(0, 0, 1) and boundary conditions (|24|1 . By setting u = 
in the equations, we obtain the conditions defining the 
quasiequilibrium state 

j • Vco = VIIo, Co = const (i), 

where the subscript "0" indicates the state of quasiequi- 
librium. Integration of the first equation leads to a con- 
dition 

RvC{){x,y,z) = Tlo{z) + J{x,y). (25) 

Here f{x,y) is an arbitrary function of longitudinal coor- 
dinates. Thus, there is a set of possible quasiequilibrium 
distributions of particles that is defined by a sum of a 
function of only transversal coordinate and a function of 
longitudinal coordinates. 

As a quasiequilibrium concentration we choose a linear 
function of the transversal coordinate 

co = a + l3z, a > 0, |/3| < a (26) 

and investigate its stability. Although gravity is unim- 
portant in the framework of our theory, in a quiescent 
fluid (in the absence of vibration) it ultimately leads to 
stratification of particles. For this reason, the choice of 
the initial distribution of particles in the form H26() is 
rather natural for a narrow layer. Further, as a reference 
value for the particle concentration we use /3/i, therefore 



B. Statement of linear stability problem 

To investigate the stability of the quasiequilibrium 
state we introduce perturbations of the velocity u, the 
pressure q, and the concentration (j) , and substitute the 
disturbed flelds u, Hq + q, cq + cj) into the problem H22(l - 
(|24|l . Wc linearize the equations near the quasiequilib- 
rium and obtain a set of equations and boundary condi- 
tions for small perturbations 

d\x 

ot 

d<j) 

— — h u • Vcq = 0, div u = 0, 
ot 

z ±1 : u = 0. 

Since the quasiequilibrium state is uniform and 
isotropic in a plane of the layer, we can restrict our stabil- 
ity analysis to the case of the two-dimensional perturba- 
tions in the form of rolls. Assuming all the perturbations 
to be proportional to exp(Ai — iky), where A = Ar -I- iXi 
is the complex growth rate and k is the real wave num- 
ber, we rewrite the problem for one variable. Eliminating 
the pressure from the equation of motion and taking into 
account the continuity equation, we obtain for the ampli- 
tude (f) of the concentration perturbation (the perturba- 
tions of the concentration and the transversal component 
of the velocity differ only by a constant factor): 

X^D^^ - R^k^4>' = XD^4>. (27) 
z = ±1 : = 0, 0' = 0, (28) 

where a differential operator D'^ — d? /dz"^ — fc^ is intro- 
duced. While obtaining H27|) . H28|) it was assumed that 
A 7^ 0. In the opposite case, A = 0, there exist obvious 
solutions (f) = const for k and (f> = (f){z) for kR^ = 0; 
the perturbations of the velocity vanish for both these 
cases. The existence of these neutrally stable solutions is 
easy to explain: although the mentioned perturbations 
deform the concentration field (|26l) , it still belongs to the 
set of quasiequilibrium solutions H25|) . 

The boundary value problem H27|) . (|28|l is a spectral 
amplitude problem, the condition of nontrivial solution 
defines the growth rate A as a function of parameters 
Rv, k. In the absence of vibration (i?„ = 0), besides 
the discussed solutions for A = 0, the problem admits 
an explicit solutionj^ eigenvalues of which belong to a 
spectrum of real negative values: all the perturbations 
decay and the system is stable. Further we are interested 
in the nontrivial case oi R^ ^ Q. In this case the problem 
(|27|l . (|28|l is no longer self-adjoint, therefore the spectrum 
of growth rates can be complex. 

Because the parameter (3 is not of the fixed sign, the vi- 
bration parameter can take on both positive and negative 
values. However, the problem is invariant with respect to 
the mutual change of sign of the transversal coordinate 
and the vibration parameter. Such a symmetry is caused 
by the absence of gravity, which allows us to restrict our 
analysis to the case i?„ > 0. 
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C. The case of small values of vibration number i?„ 

Consider the case of small intensities of vibration ac- 
tion. We look for a solution to the problem (|77|) , (^5)) in 
the form of series in small Ry 

(j) = ^o+Rv(i>i+Rl<l)2 + - . . , A = Xo+Rv\i+Rl\2 + . ■ ■ ■ 

(29) 

Substituting these expansions into Eqs. (|27|l . (|28|l . we 
find that 4>o and Ai are defined by the following problem 

XiD^^o + k^4>'o = 0, (30) 
<^o(±l) = 0, 0[,(±1) = 0. (31) 

Let us multiply Eq. H3U|I by the complex conjugate (/)g 
and integrate by parts across the layer. Taking into ac- 
count (|31|) . after straightforward transformation we rec- 
ognize that Ai — iuji, where uji is a real number. To the 
second order we obtain the problem 

XiD'^$i+k^[ = XjD^^o- X2D^$Q, (32) 

0i(±l) = O, 0'i(±l) = O. (33) 

Let us formulate the solvability condition for problem 
(Ell- We multiply Eq. ^ hy 4>*q, the solution of 
the problem adjoint to H27|l . (|28|l . and integrate across 
the layer. Taking into account boundary conditions (|31|l . 
(|33|l . we end up with 

A2 = ^ z '—^ > 0, (34) 

(|<^;,'|2 + 2fc2|^^|2+A;4|0o|2j) 

where (...) — J^^...dz. 

Thus, for any value of the wave number k the growth 
rate A2 is real and positive, and hence for any arbitrarily 
small value of the vibration parameter R^ there appears 
oscillatory instability. 



numerically. First, the spectrum of eigenvalues was ana- 
lyzed by the Galerkin method with a basis given by the 
eigenfunctions^ of the problem l|2Hl) at i?„ = 0, 

A ^ 0. Next, the behavior of the mode with the largest 
Xr was specified by the standard shooting method. '^^ It 
was found that in all the range of the parameters k and 
Ry > there exists a growth rate with the positive real 
part Xr > 0: the system turns out to be unstable for any 
intensity of vibration action. The real part of the growth 
rate as a function of the wave number (the branch with 
the maximal A,.) is plotted in Fig. ^ for different values 
of vibration parameter Ry. For any value of the vibra- 
tion parameter Ry the real part of the growth rate A^ 
has a maximum at some value /c* of the wave number. 
This critical value /c* defines the perturbation with the 
fastest growth rate, which dependence on the vibration 
parameter is given in Fig. |5] At small values of Ry the 
curve is outgoing from a point with a finite value of k^. 
With the growth of vibration intensity the characteristic 
length scale k~^ of hydrodynamic patterns developing as 
a result of instability of a quasiequilibrium state mono- 
tonically decreases. 

4.1 I ■ ■ ■ 1 
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250 500 750 1000 

FIG. 2: The critical value fc, of the wave number versus vi- 
bration parameter Rv 



D. The case of finite values of vibration number 

In the case of arbitrary values of the vibration parame- 
ter Ry the boundary value problem (|27|l . (|28|l was solved 
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FIG. 1: The real part of the growth rate Xr as a function of 
the wave number k. 



Figure 121 demonstrates the eigenfunction 0, which cor- 
responds to the breakdown of the quasiequilibrium and 
differs from the transversal velocity component by a con- 
stant factor. Under vibration action a nonuniformity 
of the quasiequilibrium distribution of particles leads to 
generation of the averaged fluid motion, which, in its 
turn, deforms the initial distribution of particles. 



IV. AVERAGED DYNAMICS OF A BUBBLY 
FLUID 

A. Basic equations 

So far, the disperse phase was assumed to be nonde- 
formable. However, the description of a bubbly fluid in 
vibration field requires a closer look. The presence of the 
disperse phase that is able to change its volume under 
the action of an external periodic field, can drastically 
influence the pulsation and the averaged dynamics of the 
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the pulsation (depending on the fast time r) and the 
averaged parts. The notation and units for the pulsation 
and averaged components of the velocities, the pressure, 
and the concentration are the same as used before, in 
Sec. m Hereafter, b is the pulsation velocity of radial 
motion. 



B. Pulsation motion 



FIG. 3: Streamlines of the averaged flow for the perturbation 
with the fastest growth rate, Rv = 10. 



system. Next, the case of a monodisperse bubbly fluid 
is studied and the criterion when compressibility of the 
disperse phase cannot be neglected is provided. Assum- 
ing the concentration ip of bubbles to be small, we formu- 
late the equations for the momentum of the liquid carrier 
phase and the bubbles Q together with Q, where for 
our particular consideration k = 0. An essential modi- 
fication that has to be taken into account here concerns 
an additional contribution to the interphase force^S 

where b is the space-averaged radial velocity of bubbles. 

Further, assuming that the size of viscous boundary 
layer i is small compared with the time-averaged radius R 
of a bubble {£ R), we can neglect the memory forces. '^^ 
Indeed, in the frames of the accepted approximation the 
leading contribution to the pulsation motion is made by 
the inertia forces [the first and the fourth terms in 0], 
whereas for the averaged motion the Stokes force and 
Fc are dominant. We suppose the density of the gaseous 
phase to be much smaller than the density of the fluid and 
neglect the left-hand side in Eq. (|2 bll . i.e., the force on the 
bubbles vanishes. The equations defining the kinematics 
readSa 



div (1 — (/j) U + div ip\Jd 



3ip 
rd 



— + (fidiv Vd = — 0, 
at Td 



(35) 
(36) 



Consider vibration of a cavity entirely filled with a bub- 
bly fluid. The principal assumptions introduced in SecHll 
are kept the same: the carrier phase is incompressible, 
the frequency of vibration is high, but subacoustic, the 
size of bubbles and the vibration amplitude are small. 
On the other hand, we impose stronger restrictions for 
the amplitude a and the frequency lo of vibration: 



niax(;y,X,Xd) 

i?2 



(39) 



where x and Xd are the thermal diffusivities of the phases. 
The second inequality implies that dissipation is unim- 
portant for the pulsation motion, therefore the pressure 
oscillation of gas obeys an adiabatic law. 

Let us now formulate the equations for the pulsation 
motion. The units for the pulsation velocity and the pres- 
sure are kept the same as earlier (see Sec. lIIB|l . as a unit 
for the velocity of the radial motion we choose aujL/R. 
The equation for the pulsation of the concentration 
splits off and can be treated separately. For this field we 
obtain 



del) ^aL 



■ div cV 



I 



R 



> 1. 



(40) 



As it follows from Eq. (|40|l . the pulsation of the concen- 
tration is small compared to the averaged field c. Note, 
that the second term in the right-hand side of Eq. H40|l 
is, generally, much smaller than the first one. However, 
under certain conditions it becomes non-negligible (see 
Sec UVDll . 

It is convenient to write down the equations for the 
remaining pulsation fields in terms of amplitudes. As- 
suming that all the fields are proportional to exp it, we 
obtain 



dpd 3 dra 

— = bpd, -TT = 0- 

at rd at 



(37) 



The pressure in the fiuid P, the pressure in a bubble Pg 
and the velocity of the radial oscillation of a bubble are 
coupled through the Rayleigh-Plesset equationi^2i2£ 



rd 



db 
di 



1 



Pn 



2cr 

rd 



P 



(38) 



where a is the surface tension. 

According to the procedure laid down in Sec. ^1 we 
introduce the time hierarchy, splitting all the fields into 



i(l - c)V = -Vg, Yd = 3V, (41) 
div(l-H2c)V = 3d^fe, q = ib{n^-l), (42) 

where V, V^, q, b are the complex amplitudes of the 
velocities of the phases, the pressure, and the velocity 
of the radial motion. Equations (|41|l . H42|l contain the 
following dimensionless parameter: 



1 



PLJ^R^ 



2(7 



(43) 



Here, J7 is the ratio of the eigenfrequency of the volume 
oscillation of a bubble to the frequency of vibration uj, 
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the parameter 7 is the adiabatic exponent, and Pg is the 
averaged pressure inside a bubble. Note the known rela- 
tion between the velocities of the phases in a vibration 
field:^^ the amplitude of oscillation of a bubble in an in- 
viscid fluid is three times as much as that of the fluid. 

Equation H42|l contains a product of two asymptotic 
parameters: large and small c. Depending on the rela- 
tion between these parameters, different solutions to the 
pulsation problem is possible. Let us introduce a renor- 
malized concentration <I>(r) = Pc (r). Assuming that this 
field is finitej^ we obtain a Hclmholtz equation for a 
potential ip of the fluid velocity and the impermeability 
boundary condition: 



3$(r) 
V„V|r 



0, 

jn- 



(44) 
(45) 



The amplitudes of the pulsation fields are expressed in 
terms of the potential ip as follows: 



-iip. 



(46) 



Equation (|44|l has a form of the stationary Schrodinger 
equation, in which the role of the scattering potential is 
played by the averaged concentration field <i>(r) - the 
waves are scattered on the nonuniformities of the bubble 
distribution. The impermeability condition (|45|l should 
be imposed on the volume- weighted velocity (1 — c)V -|- 
cVrf, which by virtue of (|41|l differs from V only by a 
small term 2cV. Note that as the frequency gets closer 
to a resonant value 0=1, the susceptibility of the system 
to vibration action gets much higher. At the point of 
the resonance, the performed analysis is no longer valid, 
because even weak dissipation becomes important and 
cannot be neglected. 

Generally, Eq. (|44|l can be solved only numerically. 
However, in the case of a small parameter e = ($) / (O^ — 
1), it is possible to obtain the leading part of the solution. 
Here, ($) is the value of <i> averaged all over the volume 
of the system. Note that the parameter e can be small 
for two different reasons: for small concentration (<i>) or 
for small vibration frequencies, O 3> 1. Assuming that 
2-axis is aligned along the direction of vibration, we can 
represent the solution in the form 

tp ^ Z - Zc + ijjl, 

where the term z — Zc corresponds to the solution of (|44|l , 
(|45|l for $(r) =0 and ~ e is defined by a boundary 
value problem 



(47) 
(48) 

Applying the solvability conditions for the problem (|47f) . 
(I48|l . we find the constant 



1 r 



(49) 



which has the meaning of the z-coordinate of the center 
of mass for the distribution of bubbles. 
Thus, for the case of small e we obtain 



V=j + 0(e), 6 = - 



02 - 1 



+ 0(e). (50) 



C. Averaged motion 

We now proceed to the averaging of the governing 
equations. The equation of the momentum for the av- 
eraged motion of the fluid takes the following form: 



\Dt J a OT 

(51) 

Here, as in Sec.m Rv — cP'liP'I? jv^ and e is the averaged 
shear rate tensor. From the averaged equation for the 
bubbles follows 

2r2 ^ 



^^Rv V • VV + i?„ 6V -f Ud - u = 0, 

where the first term is small compared to the others. 
With account of I^BJ, this equation provides a relation 
between the averaged velocities of the phases: 



Ud = u + gvV'^ 



1 Rv 



4 02 



(52) 



Here Q is the vibration parameter, which is different to 
that introduced in Sec. ITU 

To perform averaging in Eq. H51|l we take into account 
relations (|46|l and Eq. (|40|l . Retaining only the leading 
terms, we arrive at the equation 



du 

at 



u • Vu 



-vn 



(53) 



While obtaining this equation, we neglect the terms pro- 
portional to c and renormalize the pressure; the renor- 
malized variable is denoted by H. Variations of the av- 
eraged pressure, having in dimensional units the order 
pb^ ~ p{au!l)^, are assumed to be much smaller than the 
averaged pressure Pg inside a bubble. This assumption 
ensures that the averaged radius of the bubbles keeps 
constant and equals R. The absolute value of the aver- 
aged pressure becomes insignificant, only its gradient is 
of importance. 

The averaging of the equation for the bubble concen- 
tration H3()l) is straightforward and results in 



(9$ 
'dt 



+ div ^Ud = 0. 



(54) 



The right-hand side of Eq. (|54|l vanishes, because the 
phases of b and V are shifted with respect to and by 
a quarter-period. For the same reason, there is no non- 
trivial contribution to the averaged mass balance equa- 
tion, which transforms to the incompressibility condition 
for the fluid 



div u = 0. 



(55) 
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Equations together with Eqs. (gljl, igBJ defin- 

ing the pulsation potential -0, represent a set of equations 
for the averaged dynamics of a bubbly fluid. 

In perfect analogy with the case of nondeformable par- 
ticles, the averaged motion induced in a boundary layer 
is negligibly small; the principal effect is caused by the 
considered mechanism of bulk generation. Thus, at rigid 
boundaries the conventional no-slip boundary condition 
for the averaged fluid velocity should be imposed. Ac- 
cording to (|52|l . the velocity of the bubbles at the bound- 
aries does not turn to zero. Thus, the boundary condi- 
tions for the concentration of bubbles must be prescribed 
at the inflow boundaries, i.e., where n • > (n is the 
outward normal to a boundary) . 

Note that the intensity of the averaged flow of the bub- 
bly fluid is signiflcantly (by a factor of or c~^) higher 
than for the medium with nondeformable particles (see 
Sec.nj. Contrastingly to the case of nondeformable par- 
ticles, the latter circumstance enables us to consider the 
vibration parameter Ry to be finite for bubbly fluids. 
However, for low vibration frequencies, f2 ^ 1, this for- 
malism is no longer valid. To obtain any nontrivial ef- 
fects, it becomes necessary to consider the vibration pa- 
rameter Ry as asymptotically large, when the vibration 
mechanisms studied in Sec. ^become important and can- 
not be neglected any more. This situation is of special 
consideration and is studied in the next section. 



D. Low frequency approximation. Transition to a 
suspension of nondeformable particles 

The case of the low (with respect to the frequency of 
natural oscillation of a bubble) vibration frequency, when 

» 1, should be treated separately. In this particular 
situation, it is convenient to choose the value aujl/fl'^ as 
a unit for the velocity of the radial motion, the units for 
the other quantities are kept the same. As a result, the 
pulsation equations read 

z(l - c)V = -Vg, Vd = 3V, (56) 
div(l -f 2c)V = 3 nfcb, q = ib. (57) 

Here fti = fl/l is the ratio of the two asymptotic parame- 
ters. At large values of ili compressibility of the bubbles 
becomes insignificant. We note, that a similar result has 
been found out in Ref. compressibility of a bubble 
can be neglected only provided that the dimensionless 
frequency of the radial oscillation of a bubble is large 
compared to I. 

We represent the solution of the pulsation problem 
(|55|l . (|57|) in the form of series in the small concentra- 
tion of bubbles 

V = Vo + Vi + ..., q = qo + qi + ..., b^bo + h + ..., 
where Vi, gi, 5i ~ c. 



To the zero order the pulsation problem reads 

iVo = -Vgo, divVo = 0, qo = i6o, 
n-Vo|r = jn, 
which admits an obvious solution 

Vo=j, 6o = -(z-Ze). (58) 

In the same way as in Sec. IIVBI the constant Zc is ob- 
tained from the solvability condition in the next order 
and takes the form H49I) . 

As in Sec. in it is necessary to take into account a cor- 
rection Vi to the pulsation velocity. For further purposes 
only the vorticity of this field turns out to be important. 
Thus, from H5t)|l we obtain 

VxVi=Vcxj. (59) 

Equation (|40|) for the pulsations of concentration trans- 
forms into 

^ = 3^{n-Hc-yVc). (60) 

Here, in contrast to Eq. 140|l . the terms in the right- 
hand side are of the same order and therefore have to be 
retained. 

Averaging of the equations is performed in the same 
way as in Sees. Ill CI IIV CI As a result we arrive at a set 
of the averaged equations 

^ + uVu = -VH + V^u 

ot 

+R; (jj • V$ + 3^$V6^) , (61) 

divu = 0, Ud = u + R^ylV6^ (62) 

— + div$Urf==0, (63) 
ot 

where A ^ 1/i'inf), R/ = l-^Ry. 

In the limit of low vibration frequency, when ^ — + 
(ri; ^ 1), the model reduces to the case of nondeformable 
bubbles: equations H61l) - (|63() yield the same model as has 
been developed in Sec. HIl [cf. Eqs. (HSH- Indeed, 

for the case of bubbly fluid, 5 <C 1, and inviscid pulsa- 
tions, ivT > 1, we have B = 1 [see HA.2|) . (IA.3|) ] and the 
vibration force i?; jj • V'l' = Rvjj ■ Vc. 

In the opposite limiting case, A — > oo, the compress- 
ibility of bubbles becomes of crucial importance; the set 
of equations transforms to Eqs. (|52|I - H55|I . In this case, 
the solutions of the pulsation problem (|58|l and H50|l co- 
incide (if measured in the same units). 

V. EVOLUTION OF A BUBBLY FLUID IN A 
LAYER 

A. Statement of the problem 

A rather good understanding of vibration action on the 
dynamics of a bubbly fluid can be obtained from a one- 
dimensional consideration. Let us apply Eqs. H52|l - H55|) 
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to study the vibration dynamics of this medium in an 
infinite plane layer — 1 < z < 1, confined by solid walls, 
the vibration axis is transversal to the boundaries. The 
quiescent state with the uniform distribution of bubbles 
is chosen as the initial one. Symmetry of the equations, 
boundary conditions, and the initial state enables us to 
treat the problem in a half of the layer, where the func- 
tions $ and V' are even and odd, respectively. The fluid 
velocity vanishes due to the incompressibility condition 
(|55|l and the no-slip boundary conditions at the walls. 
The concentration and the potential of the pulsation ve- 
locity are defined by the following boundary value prob- 
lem: 



3$ 



-^ = 0, 



f}2_ 1 

: V = 0, $' = 0; z^l: tP' 



1, 



(64) 

(65) 
(66) 



where the prime is used to denote a derivative with re- 
spect to z. 

It is obvious, that the variation of the absolute value of 
Q results only in changing of the time scale, and without 
loss of generality, we set further Q — ±1 (the upper sign 
corresponds to SI > 1, and the lower one ~ otherwise). 
We can also renormalize the initial concentration in such 
a way that ^{z, t = 0) — 1; this leads only to rescaling of 
the frequency Q. 

We can easily make sure that there exists no quasis- 
teady state, when the fiuid and the bubbles are quiescent 
on average. Indeed, for the quasisteady state it immedi- 
ately follows from Eqs. (|65|l 



(<i>nVo^o)' = 0, 



(67) 



i.e., $0 — const/ (?/'oV'o)- As earlier, the subscript "0" 
stands for indication of the quasisteady state solution. 
On the other hand, Eq. 1)64(1 then transforms to 



3 const 
- 1) 



0. 



(68) 



Integrating this equation and taking into account that ipQ 
is even, we conclude that const — and hence $0 = 0. 
Note, that this conclusion becomes invalid at the points 
where ipQ'4'o — 0. 



B. Initial stage of evolution 

We now turn to investigation of the initial stage of 
evolution, when an analytical solution can be obtained. 
We assume, that $ = 1 + ^i{z,t), where |$i| < 1. 

Neglecting the small perturbation $i in the 
Schrodinger Eq. (|54|l . we have 



^0 = 0. 



(69) 



First, consider the case fl > 1, when the vibration 
frequency is lower than that of the natural oscillation of 
a bubble. We introduce 



a 



- 1 



(70) 



and note that this expression for the wave number a is 
in agreement with the known disperse relation for waves 
in a bubbly medium (see, for example, Ref. 4Qi4ll) . 

The solution to Eq. (|69|l with the impermeability con- 
dition at the boundaries is as follows 



sm az 
a cos a 



(71) 



The resonant frequencies correspond to the divergency of 
()71|l and are defined by the expression 



ni = 1 



12 



7r2(2n+ 1)2 



0,1,, 



(72) 



Note, that at these frequencies it is necessary to account 
for dissipation in the pulsation equations. 

Taking into account the result H71|l . we obtain the av- 
eraged velocity of the bubbles and the correction to the 
concentration 



Ud 



sin 2az 



a cos a 



$1 = 



2 cos 2az 



■ t -\- const. 



(73) 



We see that the most rapid decrease of the concen- 
tration of bubbles occurs in the vicinities of the points 
Zn = -Kin/a^ n ^ Z - the bubbles tend to leave the 
nodes of the pulsation pressure; and vice versa, the 
concentration grows in the antinodes of the pressure. 
This phenomenon is referred to as the primary Bjerknes 
effect However, we stress, that in our case, the bub- 
bles are not just advected by the external nonuniform 
pulsation field (as in the conventional Bjerknes effect). 
The presence of the bubbles causes this nonuniformity. 

Let us proceed to the opposite range of frequencies, 
SI < 1. In this case, the frequency of external action is 
beyond the passband for the bubbly fluid, i.e., the waves 
decay deep into the layer: 



sinh ( 



sinh2az 

u ' " ^ 1 7^' '^d = -2 — , (74) 

a cosh a 1 - S2^ cosh a 

Further, for the perturbation of the concentration we find 



$1 



2 cosh2Q!Z 



cosh 



t + const. 



(75) 



Thus, the bubbles migrate away from the boundaries 
of the layer to its center. As it could be expected, at a 
frequency larger than the resonant frequency of a bubble, 
the bubbles accumulate in the node of the pressure, at 
the center of the layer. 
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C. Finite time evolution. Numerical results 

At finite times the set of equations H52l) - (|55|) . describ- 
ing one-dimensional dynamics of a bubbly fluid in the 
layer, was studied numerically. The equation for transfer 
of the bubbles was integrated by means of the method of 
characteristics. The equation for the characteristics has 
an obvious form 



dz 
It 



Ud{z,t). (76) 
At the characteristics the following relation is fulfilled 



(77) 



The fields of the concentration <i> and potential tp are 
defined in the nodes of the uniformly spaced grid; we 
used up to 2000 grid nodes. In the case when Mrf(l) < 
0, there appears a front of the concentration z — Zs{t) 
such that $ ^ only in the domain with z < Zg. For 
this particular situation, the grid nodes were interposed 
only in the domain occupied by the bubbly fiuid, i.e., for 
z < Zs(t); in the opposite part of the layer, the equations 
admit an obvious solution: ^ — 0,ip — ipo + z. At 
the point z — Zg the functions ip, "0' £^re continuous, so 
that the value of the constant ipo is insignificant. In the 
opposite case, Ud{l) > 0, a part of bubbles settles on the 
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FIG. 4: The maximal ^max and the averaged (<E>) values of 
concentration as functions of time (a); the concentration pro- 
files at different times, Q — =3 (b). 
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FIG. 5: The maximal ^max and the averaged (4>) values of 
concentration as functions of time (a); the concentration pro- 
files at different times, Q — 1,Q^ = 1.25 (b). 

boundaries: the averaged concentration ($) = ^ dz 
decreases with time. 

Results of the numerical simulation are presented in 
Figs.EEl If > fio = 1-489 [for the uniform distribution 
of the bubbles, see relation JT^J], the bubbles gradually 
leave the bulk of the layer, settling on the boundaries, 
Fig. 0] At first, the maximal value $m of the concen- 
tration grows and after a while starts to decrease; the 
averaged concentration (<f>) monotonically decreases with 
time. 

In the case of fii < < JIq, the bubbles accumulate in 
a certain domain inside the layer. The concentration of 
the bubbles rises abruptly near the point z^ and reaches 
infinite values for a finite time t — to - there develops 
a peaking regime. The maximal concentration grows as 
(t — to)~^, i.e., a bubble screen develops. The potential 
itself remains finite, while its second derivative tends to 
infinity near the point Zs- This scenario is demonstrated 
for il^ = 1.25 in Fig. [S] The averaged concentration of 
the bubbles decreases with time. As the eigenfrequency 
O gets closer to 1, the number of bubble screens increases. 

In the case < 1, the bubbles accumulate in the center 
of the layer - a node of the pulsation pressure. Evolu- 
tion of the maximal concentration, the coordinate z^ of a 
front, i.e., the interface bubbly fluid-pure fluid, and the 
profiles of the concentration are depicted in Fig. As 
was mentioned above, the bubbles located near the front 
possess the maximal velocities [cf. (|7^ ]. consequently 
the concentration grows very fast at the front and does 
not practically vary in the stagnant core in the center of 
the layer. 
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FIG. 6: The maximal concentration ^max and the coordinate 
of the front Zs as functions of time (a); the concentration 
profiles at different times, Q = —1, = 0.5 (b). 



VI. CONCLUSIONS 

The averaged dynamics of various two-phase systems 
in a high-frequency vibration field has been theoretically 
studied. The continuum approach was applied to de- 
scribe such systems as solid particle suspensions, emul- 
sions, bubbly fluids, when the volume concentration of 
the disperse phase is small and gravity is insignificant. 
The dynamics of a monodisperse system was considered 
by means of the averaging method, when the fast pulsa- 
tion and slow averaged motion can be treated separately. 
The vibration dynamics of suspensions of both nonde- 
formable and deformable particles, when the compress- 
ibility of the disperse phase becomes important, has been 
investigated. 

An averaged model for nondeformable particles has 
been obtained. It is shown, that accounting for the Stokes 
drag and unsteady forces such as the Basset history term, 
the second memory force and the added mass force in the 
pulsation motion, is of crucial importance for correct de- 
scription of the averaged dynamics. The averaged equa- 
tions have been obtained and simplified in the framework 
of the single-fluid approximation. It is proven that the 
particles can be treated as frozen into the averaged fluid 
fiow. As a result the averaged dynamics of a suspension 
is described by the equations of momentum and mass 
conservation for the carrier fluid and a transfer equation 
for the concentration of particles. The action of vibra- 
tion results in the appearance of a vibration force in the 
equation of the fluid motion. This force is nonvanishing 



only for the nonuniform distribution of particles, its di- 
rection coincides with the vibration axis. We note, that 
the developed model is expected to be worth applying to 
clarify recent experiments on blood flow resistance under 
vibration action. 

The developed model has been applied to study the be- 
havior of a two-phase medium in an infinite plane layer 
subjected to transversal vibrations. It is demonstrated 
that there is a set of possible nonuniform quasiequilib- 
rium distributions of particles, when the averaged fiow 
vanishes but the pulsation velocities do not. The stabil- 
ity of a linear distribution of particles has been investi- 
gated. Analytical and numerical analysis manifests that 
the quasiequilibrium state is unstable for any intensities 
of vibration. A nonuniformity of particle distributions 
causes the vibration force that generates the averaged 
motion in fiuid. Because of this fiow, the initial distri- 
bution of particles is deformed and the state becomes 
unstable. 

As a system with deformable particles, the case of bub- 
bly fiuids in vibration field has been analyzed separately. 
An averaged model accounting for the compressibility of 
the bubbles has been developed. The pulsations were as- 
sumed to be inviscid, and a set of averaged equations has 
been obtained in the single-fiuid approximation. The av- 
eraged velocity of bubbles differs from the fluid velocity, 
this difference is caused by vibration and proportional to 
its intensity. It is shown that in contrast to the case of 
nondeformable particles, the impact of vibration on the 
system with deformable particles is signiflcantly stronger. 
Even for uniform distribution of particles the vibration 
force is nonzero. 

An intermediate case of low vibration frequencies, 
when the ratio f2 of eigenfrequency of radial bubble oscil- 
lations to the the frequency of vibration is high, has been 
considered. A criterion when the compressibility of bub- 
bles can be neglected has been figured out: Q ::§> L/R, 
where L is the length scale of the flow and R is the time- 
averaged radius of the bubble. In this case the inter- 
mediate model reduces to the model for nondeformable 
particles. In the opposite limit, 51 <C L/R, the inter- 
mediate model transforms to the discussed model where 
compressibility of bubbles becomes of crucial importance. 

The dynamics of bubbly fiuid in an infinite plane layer 
under the action of transversal vibration has been ana- 
lyzed. The quiescent state with the uniform distribution 
of bubbles is chosen as the initial one. It turns out that 
for > 1.489 the bubbles migrate to the boundaries of 
the layer; for 1 < < 1.489 the bubble screens appear; 
for < 1 bubbles accumulate in the center of the layer. 
We have demonstrated that the behavior of bubbly fiuid 
is analogous to the primary Bjerknes effect: for Q > 1 
the bubbles leave the nodes and accumulate in the antin- 
odes of the pressure wave, while for Q < 1 the bubbles 
migrate to the nodes. However, in our case, the bubbles 
are not only advected by the external nonuniform field, 
as in the conventional Bjerknes effect, but also cause this 
nonuniformity. 
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To demonstrate the dependence in more general situ- 
ations, we tabulate (|21|l as a function of K for different 
values of 5 and n. Typical results are presented in Fig.[7| 
As it can be shown, the function B{K) is positive for 
all K] it is either monotonic or has a maximum. In the 
limits of A' ^ 1 and A ^ 1 the results are independent 
of K and agree with the expressions ljA.l|l and (|A.2|I , re- 
spectively. As it can be expected, the effect of different 
values of k occurs at intermediate values of K . 



APPENDIX 

Here we analyze relation (|21|l as a function of the pa- 
rameters (5, K and K. First, we give the expressions for 
the limiting vibration regimes of viscous (A <C 1) and 
inviscid (A ^ 1) pulsation. Next, we discuss the de- 
pendence in the whole range of values of K and outline 
features for some typical media. 

In the limit of viscous pulsation, A ^ 1 , relation H21() 
reduces to the following: 



So 



(A.l) 



Note, in this limit, it is not enough to take into account 
only the Stokes force, which would result in vanishing 
B and therefore in no vibration force. Physically, this 
means that the particles get frozen not only in the av- 
eraged flow, but in the pulsation fluid flow as well, and 
therefore no averaged effects are possible. To account for 
any nontrivial dynamics one has to retain small (with re- 
spect to the Stokes drag) terms due to memory forces, 
which ensures relative pulsation motion of phases. An- 
other interesting point is that this result does not depend 
on Hi. 

In the opposite limit of inviscid pulsation, A 3> 1, the 
leading contribution to B is made by the added mass 
force. The result is sensitive solely to the relative density 
of phases: 



^- 1 
1 + 25 



(A.2) 



A typical case of a bubbly fluid (k ^ 1, (5 ^ 1) is 
governed by 



Bh^ 



2A3(18\/2 + 9A 3^/2A2 + A^) 



(18V2 + 18A - A3)2 + A2(18 + 3V2A + K^Y ' 

(A.S) 

The result is caused by the Stokes drag, the second mem- 
ory force, and the added mass force. 

For the case of solid particles (k ^ oo), the dependence 
(|21|l transforms to: 



B. = 



((5- l)2^/2A3(9 + ^/2A) 
81(\/2 + A)2 -f A2[9 + V2(l + 25)Y ' 



(A.4) 



where the Stokes drag, the Basset and the added mass 
forces make finite contribution. 
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FIG. 7: The dependence of B against the parameter K for 
5 = (a) and (5 = 4 (b). 



Figure [7{a) represents the case of bubbly fluids {5 <^ 
1); the curve with k = corresponds to relation (|A.3p . 
The curves go out from the point i? = at A = 0, start 
to grow at small A as A'^/(9\/2) [see relation (jA.ip ]. and 
monotonically approach the value Boo = 1, as follows 
from ljA.2|l . For S < 1/2, qualitative form of the curves 
keep the same as in Fig. [TJa), the value of Boo changes 
according to ljA.2p . 

For 6 > 1/2 the dependence B{K) is no longer mono- 
tonic: there appears a maximum Bm at some value A,„ 
[see Fig. [3b)]. The curve with k ^ 1 corresponds to 
solid particles and is described by formula HA.4|I : such a 
curve in Fig. E^b) is plotted for a suspension "sand in 
water" ((5 = 4, k 3> 1). We note, that for higher rel- 
ative densities 6 the dependence B{K) does not change 
qualitatively. In the limiting case of heavy particles sus- 
pended in a gaseous medium (solid particle suspension in 
air, aerosol), 6^1, the maximum value Bm grows and 
shifts to lower values of A; the asymptotic behavior reads 
Bm oc ^/S, Km oc 1/ ^/S. At high values of A parameter 
B gradually approaches the limiting value Boo = 1/4. 
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